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QUANTILE CLOCKS 

By Lancelot F. James-*^ and Zhiyuan Zhang 
Hong Kong University of Science and Technology 

Quantile clocks are defined as convolutions of subordinators L, 
with quantile functions of positive random variables. We show that 
quantile clocks can be chosen to be strictly increasing and continuous 
and discuss their practical modeling advantages as business activity 
times in models for asset prices. We show that the marginal dis- 
tributions of a quantile clock, at each fixed time, equate with the 
marginal distribution of a single subordinator. Moreover, we show 
that there are many quantile clocks where one can specify L, such 
that their marginal distributions have a desired law in the class of 
generalized s-self decomposable distributions, and in particular the 
class of self-decomposable distributions. The development of these 
results involves elements of distribution theory for specific classes 
of infinitely divisible random variables and also decompositions of a 
gamma subordinator, that is of independent interest. As applications, 
we construct many price models that have continuous trajectories, 
exhibit volatility clustering and have marginal distributions that are 
equivalent to those of quite general exponential Levy price models. In 
particular, we provide explicit details for continuous processes whose 
marginals equate with the popular VG, CGMY and NIG price mod- 
els. We also show how to perfectly sample the marginal distributions 
of more general classes of convoluted subordinators when L is in a 
sub-class of generalized gamma convolutions, which is relevant for 
pricing of European style options. 

1. Introduction. Let Qr{u) = uii{t\ Fji{t) > u},0 < u < 1 denote the 
quantile function of a nonnegative continuous random variable R with strictly 
increasing cumulative distribution function (c.d.f.) Fr, and finite first mo- 
ment E[i?]. In this paper, we introduce and describe detailed distributional 
properties of a class of random time changes Tr := {Tji{t),t > 0), which we 
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call quantile clocks. These processes are defined as, 

(1.1) TR{t)= f QR{{l-s/t)+)L{ds), t>0, 

Jo 

where L is a subordinator, and (a)+ := max(0, a). While applicable in many 
settings, we follow the framework in [4] and discuss the modeling advantages 
of quantile clocks as business activity times in time changed models for asset 
prices. 

The quantile clocks, may be written as special cases of convoluted subordi- 
nators, which are processes described in Bender and Marquardt [4]. That is 
to say, processes defined as T{t) = k{t, s)L{ds),t > 0, for k(t,s) a known 
kernel. The authors [4], Proposition 1, provide mild conditions on k{t,s) and 
L such that a process T := {T{t),t > 0) has almost surely strictly continu- 
ous and increasing sample paths. In terms of applications, [4] argue that one 
can use W{T(t)), where is a brownian motion with drift, as time changed 
models for the log price of assets that possess continuous trajectories, where 
T{t) is now interpreted as business activity time. Furthermore, such models 
can correct deficiencies in Black-Scholes type price models. In particular, it 
is known that (i) the log returns of asset prices have nonnormal distributions, 
and often exhibit semi-heavy or heavier tail behavior, (ii) the volatility or 
variance is dependent on time, (iii) asset prices exhibit volatility clustering 
or persistence. Reference [4] also describe a general formula for European 
style option prices that depend on the marginal distribution of T(r) for 
some fixed time to maturity r > 0. For other applications of processes repre- 
sentable as convoluted subordinators, see, for instance, [18, 30, 31, 37] and 
references therein. 

In the literature, exponential Levy price models, defined as e~-^^*^ for a 
Levy process x on M, have been quite successful in terms of their ability to 
capture some of the stylistic features of asset prices (i) and (ii) listed above. 
In addition, there are many choices of x where one can easily calibrate 
pricing models to the options market, capturing volatility smiles and skews, 
via Monte Carlo methods or perhaps more generally by the fast Fourier 
transform (FFT) methods outlined in Carr and Madan [10]. Many Levy 
processes X; can be expressed as W{C{t)) for some subordinator (. However, 
the precise ^ that is associated with a x is not always known explicitly, and 
X is often modeled via its Levy density. Arguably, the most popular models 
of this type include the variance gamma process (VG) by [28], where C is 
a gamma subordinator, the Carr-Geman~Madan-Yor (CGMY) process [7], 
where ^ has only recently been identified by Madan and Yor [29], and the 
normal inverse Gaussian (NIG) process [2], where C is an inverse Gaussian 
process. The popularity of these models is due to their relative simplicity 
and distributional flexibility. 
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However, since x has independent increments, exponential Levy processes 
are unable to capture effects due to volatility clustering. One approach dis- 
cussed in [8], which is related to the price models in [3], is to further time 
change x by a stochastic volatility process of the form T{t) = jQv{s)ds, 
where v{s) represents the instantaneous volatility either following a mean 
reverting Cox-Ingersoll-Ross (CIR) process or a non-Gaussian Ornstein- 
Uhlenbeck (BNS-OU) model of Barndorff-Nielsen and Shephard [3], speci- 
fied by the dynamics 

(1.2) dv{t) = -Xv{t)dt + 'd{Xdt), 

where "d is a subordinator we shall call an OU-BDLP. The BNS-OU model 
v{t), possesses jumps, and has a stationary distribution with laws equating to 
the class of laws of self-decomposable random variables that remarkably one 
can choose based on a prescribed choice of "i?. This latter fact is important 
for our exposition. In contrast, the integrated volatility T{t) = v{s) ds is 
continuous and has nontypical marginals laws (obviously depending on t). 
In the case of the CIR process, v{t) is a diffusion having a transition density 
following a noncentral chi-squared distribution. We will not consider models 
of this type. 

The authors [4] demonstrate that their approach, involving convoluted 
subordinators, can be viewed as viable variations of the idea in [8]. Further- 
more, their work essentially contains the popular model of [3]. However, as 
noted by the authors, there are practical issues arising for instance in the 
pricing of options that relate to the marginal distributions of T(r) at matu- 
rity times r. While T(r) are infinitely divisible, their marginal distributions 
and also characteristic functions depend on k{T,s) and r, in a nontrivial 
way. Hence, leading in general to nonfamiliar distributions for T{t). Re- 
lated to this point, is a classical problem where in general it is not clear 
how to exactly sample infinitely divisible random variables, even in the case 
of a subordinators ^(t) for each fixed r. Some notable exceptions for ( are 
gamma, positive stable and inverse Gaussian processes whose marginals for 
each fixed time point are gamma, positive stable and inverse Gaussian ran- 
dom variables, and hence are easily sampled. More generally, one can resort 
to sampling methods based on truncation of infinite series representations, 
but these do not yield exact samples and it is not always clear how to control 
the level of accuracy. Reference [4] do point out that if T(r) has an ana- 
lytically tractable characteristic function then one can apply the popular 
fast Fourier transform (FFT) techniques in Carr and Madan [10] to obtain 
explicit option prices. They also provide supporting results for some choices 
of k and L. However, due to the generality of k{t,s), it appears difficult to 
apply the (FFT) for option price formula depending on general T(t). These 
points do not refiect a deficiency in the approach of [4] but rather that the 
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class of convoluted subordinators is quite general. The task then becomes 
how to choose kernels k and subordinators L that are convenient in terms 
of implementation as well as having general modeling flexibility. 

1.1. Contributions and outline. In Section 3, we show that quantile clocks, 
which arise by setting k{t,s) = Qr{{1 — s/t)+), can be chosen to have con- 
tinuous trajectories and in general have marginal distributions that, for each 

fixed time t, equate to a single subordinator ^. That is Tji{t) = (^{t) for fixed 
t. We also highlight a very tractable example related to [20, 22]. Of course, 
in general, the law C, that is, the marginal laws of Tr, depends on {R,L), 
and hence the deterministic quantile function Qji. However, in Section 4, 
we show that there are many quantile clocks whose marginal distributions 
can be chosen such that they do not obviously depend on Q/j. In fact for 
these given Qr, we show that one can choose L, and a random variable Y, 
such that the marginals of Tr have specific laws in the Jurek's [17, 23-25] Us 
class of generalized s-selfdecomposbale laws, for 6 > 0. These classes contain 
the important class C of self-decomposable distributions on M+. See [9] for 
the relevance of self-decomposable Levy processes in financial modeling. 

This ability to choose specific (familiar) marginal laws for price processes, 
while allowing for quite varied path properties induced by different quan- 
tiles Qn, gives modelers a great deal of flexibility. It is also reminiscent of 
how one might choose a BNS-OU model v to have a specific stationary dis- 
tribution, that is, v{0) = v{t) for all fixed t, in C based on the OU-BDLP 
1? appearing in (1.2). However, recall that v has jumps and the law of Tr 
obviously must depend on t. The precise methods we use to establish these 
results, and identify Qr, Y and L, are given in Sections 4, 5 and 6, and 
should also be of general interest to experts in Levy processes. In Section 7, 

we exploit the fact that Tji{t) = C(t) for each t, and we show that compo- 
sitions (or time changes) involving quantiles clocks behave marginally like 
subordinators. In Section 8, we show that as consequences of our results, 
that we are able to identify price processes whose marginal behavior co- 
incides with those of exponential Levy price processes. In particular, we 
identify explicitly many processes whose marginal distributions are equiv- 
alent to VG, CGMY and NIG price processes, but whose trajectories are 
continuous and otherwise quite varied, and additionally exhibit volatility 
clustering. We also identify models possessing jumps that otherwise have 
the properties mentioned above. In Section 9, we show how one can use our 
results for quantile clocks to specify laws for the convoluted subordinator 
referred short memory kernel in [4]. 

While quantile clocks are our main focus, in Section 2, we also describe re- 
sults that apply to the practical implementation of log price models W{T{t)), 
considered in [4], where T(t) is based on a general kernel k{t, s). In particlar, 
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if L is chosen to have laws in the class of generalized gamma convolutions 
with finite Thorin measure, see [6, 19, 22], call this class then the ran- 
dom variable T(r) can be exactly sampled in many instances. This is based 
on a very recent work of Devroye and James [15] where a double coupling 
from the past (Double CFTP) perfect sampling routine is devised, and also 
results described in James [19]. Furthermore, for this choice of L, by using a 
deterministic time-change we obtain a simplified version of the option price 
formulae given in [4]. 

1.2. Preliminaries. We now present some concepts and notation we shall 
use throughout. First, for fixed positive numbers (a, 6), let 7a denote a 
gamma(a) random variable with shape parameter a and scale 1, let f3{a,b) 
denote a beta (a, 6) random variable. Furthermore, U will always denote 

a Uniform[0, 1] variable, and recall that for any 6 > 0, U^^^ = /3s^i. is 
a Bernoulli random variable with success probability p. In addition for a 
generic random variable Y, Y' will denote a variable equivalent in distri- 
bution but otherwise independent. {N{s) : s > 0) will denote a homogeneous 
Poisson process with intensity E[A^(s)] = s. For a (nonrandom) function 
g{x), g'{x) and g"{x) with denote its first and second derivatives. 

Formally, recall that a subordinator <^ = t > 0), is an increasing pro- 
cess with right continuous paths and stationary independent increments, 
whose law is specified by its Laplace transform for some w > 

(1.3) E[e~'^^(*)] =e-*^^('"\ 

where for some c > 0, cj > 

'ip^{u) = cuj+ / (l-e-'^")A^(ds) 

is finite and is called the Laplace exponent of A^((is) is its Levy measure, 
Pc_{s) = h.(^{ds)/ds is the Levy density. We will work with the case where 
c = 0. 

It follows that the laws of C, can be specified by any of these quantities. 
Q is said to be of infinite activity if A^(oo) = oo and otherwise of finite 
activity. In the latter case, C, corresponds to compound Poisson process whose 
jumps have a common probability density/mass function proportional to p^. 
Throughout we shall reserve the notation L, Z for generic subordinators, 
and corresponding random variables, and the notation -d for the OU-BDLP. 
As is well known, for each fixed t, <^(t) is a random variable in the class J of 
infinitely divisible random variables (taking values in M+). We now describe 
the characteristics of some important subclasses of J, say C,B,Q and 
satisfying g+ CQ C C dj and g+ dj . 

We say that a random variable C(l) is in the class C of self-decomposable 
variables if P(^{s) = s~^h{s), with h decreasing. We also note that from Jurek 
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and Vervaat [26] that, with respect to the OU process in (1.2), there is the 
relationship 

C{l)£v{0) = v{t)= / e-M(ds). 
Jo 

Note that we will say that C is a subordinator in C to mean that it is a 
subordinator whose Levy density corresponds to that of a variable in C, 
of course C(^) is in ^ for each fixed t. Similar statements will apply for 
other classes. We say that ("(1) is a variable in Bondesson's [6], Section 9, 
B class, or the class of generalized convolutions of mixtures of exponential 
distributions (GCMED), if the Levy density is completely monotone, that 
is, pc^{s) = e~^y fj,{dy), for some nonnegative measure jj,. 

We now describe the classes Q and Q^. C(l) is a variable that is a general- 
ized gamma convolutions (GGC), see [6], if it is in the class G, characterized 
by 

roo /"OO 

p^{s)=s-^ e^'^uidy) and tP^{uj) = \og{l + u /y)v{dy) 
Jo Jo 

for some sigma-finite measure u, formally known as a Thorin measure. We 
say that C is a GGC(i/) subordinator. 

A C(l) variable is in the class G^, if it satisfies 

(1.4) p^{s) = 9s-^E[e-'/^] and ip ^ (lo) = eE[\og{l + loR)] 

for some > and some random variable R satisfying E[log(l +a;-R)] < oo. In 
this case, we say C(l) is a GGC(^, R) variable. Moreover, ({t) is a GGC{9t, R) 
variable for each fixed t, and C is referred to as a GGC{0,R) subordinator. 

We now highlight some important properties of GGC(^,i?) random vari- 
ables, and subordinators that for instance allow them to be exactly sampled 
by the methods in [15]. These facts can found in [19] as well as [20, 22], and 
depend heavily on the results for Dirichlet means in [13] . Letting Zg denote 
a GGC(^,-R) subordinator, it follows that ZQ{t) = Zgt{l). Importantly, there 
is the representation, for any k > ^ > 0, Zq{1) = ^oMq = "f^Mn where 

(1.5) Me = Pe,iMe + (1 - Pe,i)R 
and 

(1.6) MAPK,iM, + {l-^.,i)Rip 

for p = 9 / n. That is, a GGC(^, R) random variable variable is a GGC(«;, R^p) 
variable. In particular, if < 0t < 1, then Zgt{l) = Zg{t) = 71M1, where 



(1.7) 



Ml = UMi + (1 - U)R^p 
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for p = 6t, and it follows from the work of Cifarelli and Regazzini [13] that 
Ml has density of the form 

^ sinfvrF^^ (x))e-P*«(^) for x > 

with 

=E[log|x-i?|I(^^^)]. 

Thus, as pointed out in [19], if one one can evaluate ^'_r(x) in a suitable 
fashion, then one can exactly sample any variable Zg{t) for every fixed 
< t < 1/9, by, for instance, rejection sampling. Since any number s > 0, 
can be set to s = nt, for some integer n and < t < 1/6, it follows that 

Zq{s) = Zo^tn) can be exactly sampled by at most exactly sampling n copies 

of the random variable ZQ{t) = 71 Mi. We note that in general Mg for 9 > 0, 
does not have a simple expression for its density. So the exact sampling 
method suggested above relies solely on the ability to sample the variable 
in Ml in (1.7), for each p. This is possible provided that ^r{x) is ana- 
lytically tractable. However, since R can be quite arbitrary this will not 
always be true. Fortunately, there is the recent Double CFTP perfect sam- 
pling method by [15] that can be used to exactly sample any of the variables 
satisfying (1.5), (1.6) or (1.7). This procedure applies provided that i? is a 
bounded variable and one has a method to sample R, but otherwise does 
not require any potentially complicated calculations. Hence, any GGC{9,R) 
variable, with R bounded, can be exactly sampled by drawing an indepen- 
dent gamma variable and applying the Double CFTP. Details may be found 
in [15] , however we shall sketch out the details for a subclass of the variables 
T(t) in the next section. 

Remark 1.1. Letting Qr{u) denote a quantile function of R, variables 
Mq, satisfying (1.5), are called Dirichlet means since they can always be 
represented as 

Me= [' QR{u)Do^e{du\Fu)= H yDe{dy\Fti), 
Jo Jo 

where 

CXD 00 fc — 1 

Dg{y\Fn) =Y.Pk\R,<y) = 11(1 - y;)V,<,) 

k=l k=l j=l 

is a Dirichlet process with (P^) a sequence of probabilities having a Poisson 

Dirichlet law with parameter 9, see [16, 33]. That is, for each k, 'ygPk = Jk, 
where (Jk) are the ranked jumps of a gamma(0) subordinator. (Vfc) are i.i.d. 
Beta(l,^) random variables, and (Rk) are i.i.d. Fr. See [22, 27] for more 
details. 
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Remark 1.2. For < a < 1, positive stable subordinators Sa{t), where 
5q,(1) := Sa, with ^Sai^) = w"^, and corresponding processes Sa{t), with 
'i/'g; (oj) = (1 + a;)" — 1, as well as their scaled variations, are in Q but not G^. 
Naturally a gamma(0) subordinator, say {'ye{t)]t > 0), is in However, 
Sa,Sa and 7^, constitute a family of (generalized gamma) subordinators 
with Levy density 

for < a < 1 and 6 > 0, see [33], Proposition 21. Additionally, heavy tailed 
variables such as Linnik variables of the form Sai'^et) = iK^^a are in Q^. As 
well as their exponentially tilted counterparts Sa (jetp) , for some < p < 1 . 
See [20]. 

2. Convoluted subordinators. We now give the formal specifications for 
convoluted subordinators as defined in Bender and Marquardt [4]. Through- 
out the rest of the paper, let L denote an infinite activity subdordinator. 
That is the Levy measure, Ai(oo) = 00. In order that the convoluted sub- 
ordinator 

(2.1) T(t)= [ k{t,s)dL{s) 

Jo 

has strictly continuous and increasing trajectories, k is chosen to satisfy the 
following regularity conditions: 

(a) for fixed t £ [0,oo), the mapping s 1— >• k{t,s) is integrable, 

(b) for fixed s € [0,oo), the mapping 1 1— )• k(t,s) is continuous and in- 
creasing and there is an e > such that 1 1— t- k{t, s) is strictly increasing on 
[s,s + e], 

(c) k{t, s) = whenever s > t>0. 

The authors also derive a weighted Black-Scholes pricing formula for Eu- 
ropean style options as follows. Let 

(2.2) Wf,{t) = W{t) + nt 

denote a standard Brownian motion with drift parameter fi, that is W{t) 
is a standard Brownian motion. Recall that for geometric Brownian motion 
the price process under the risk neutral measure is given by 

(2.3) Sit) = 5(0) exp{rt + W_,/2{aH)}, 

where TV_i/2 is defined by (2.2) with // = —1/2. 

Setting Sr = S{t), the quantity (Sr — K)+ is the the payoff function of 
a European call option with strike K > and maturity r, and r > is the 
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risk-free interest rate. Then the Black-Scholes formula for the price at time 
0, say B{a,K,T), is given by 

(2.4) B{a, K,t) = e~^-E[{Sr - K)+] = SoHdiia)) - Ke-'-^<^id2{a)), 

where ^>(x) is the standard normal distribution function 

log(go/K) + (r + aV2)r ,,,,,,, r 
di[a) = — and d2{(T) = di[a) — ay/r. 

See Schoutens [35] for this notation. In [4], a price model under risk neutral 
dynamics is specified as 

(2.5) 5(t) = 5oexp{rt + Ty_i/2(a2T(t))}, 

where now S is the asset price, and T(t) is a convoluted subordinator. They 
obtain the following pricing formula. 

Theorem 2.1 (Bender and Marquardt [4], Theorem 4). For the 'price 
model (2.5), with Sr = S{t). Let {Sr — K)+ be the payoff function of a 
European call option with strike K € and maturity r. Then the initial 
fair price of {Sr — K)^, is given by 

(2.6) e-'^E[{Sr - K)+]=E[Biay^T{Tj/^,K,T)], 

where for positive y, B(y, K, r) is the Black-Scholes price given in ( 2.4 ) with 
y in place of a . Furthermore, Sq is considered fixed. 

As noted in [4], and discussed in the Introduction, the problem with the 
above result is that it is in general difficult to handle the exact law of T{t). 
However, the authors do point out that if T{t) possesses an analytically 
tractable characteristic function then it is possible to use fast Fourier trans- 
form (FFT) methods. They give some special examples where this might be 
possible, but in general this is not straightforward. This is clear since the 
Laplace exponent of T(r) can be expressed as 

(2.7) lPT{r){u:)=TE[^PL{uJk{T,UT))], 

where k and L are quite general and the expression otherwise depends on r 
in a nontrivial way. 

We believe Theorem 2.1 does have quite a bit of utility provided that one 
can have more control over the choice of marginal laws exhibited by T(t), 
for each fixed r. Next, we show that by choosing L to be in one can (in 
a practical sense) use Theorem 2.1 for many kernels k. Even those that do 
not admit nice characteristic functions. 
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Remark 2.1. We note that the martingale argument used in [4] is dif- 
ferent than that used for standard time changed models. The filtration used 
by [4] preserves the martingale property for a larger class of models includ- 
ing, of course, time changes by a simple subordinator. However, the usual 
filtrations used for simple subordinators may not preserve the martingale 
property for all convoluted subordinators. 

2.1. A general result for L in Q^. As we just mentioned, we now look at 
the choice where L is a GGC{9,Y) subordinator where Y is some random 
variable. In terms of modeling for general T{t), we can view as a time 
parameter that can be manipulated for practical convenience. This is partly 
because the variable Y can have unknown parameters that can be used for 
calibration. 

Theorem 2.2. Let W^{t) denote a Brownian motion with drift as de- 
fined in (2.2). Let 

(2.8) T{t)= [\{t,y)Le{dy) 

Jo 

denote a convoluted subordinator where Lg is a GGC(^,y) subordinator. 
For each fixed t, define a random variable Rt = k{t,Ut). Then the process 
{Wfj^{T{t)) :t>0) is almost surely continuous and has the following distri- 
butional properties: 

(i) for each fixed t, T{t) is a GGC{6t, RtY) random variable satisfying 
T{t) =^gtMet, where 

(2.9) Met = (3et,iMet + (1 - (3et,i)RtY, 

(ii) ifO<et = p<l, then T{t) = 7iMi,t , where 

(2.10) M,,t = UM^,t + {l-U)RtYCp, 

(iii) the density of the Mi^t is given by 

^sin(7rFt(x))e-P**(^) for x > 0, 
vr 

where Ft = Fji^y^p is the c.d.f. of the variable RtY^p and 
^t{x)=n^og\x-RtY\\R^Y^,-{\. 

Proof. Following (2.7) and (1.4), it is easy to see that the Laplace 
exponent of T[t) is given by 

V'T{t)(w) = tm^Le{^k{t, Ut))] = et¥.[\og{l + ujRtY)]. 
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The results (i), (ii) and (iii) then follow from the material we discussed at 
the end of Section 1.2. □ 

We now state a result for European style options, which is immediate 
from Theorems 2.1 and 2.2. 

Proposition 2.1. For the price model (2.5), let T{t) he the process 
specified by (2.8) and otherwise consider the setup in Theorem 2.1. Let {Sr — 
K)^ be the payoff function of a European call option with strike K G M_|_ and 
maturity r. Then the initial fair price of {Sr — K)^, is now given by 



(2.11) e-'-E[{Sr-K)+]=E[B{a^jgrMe,r/r,K,T)]. 

The pricing formula in (2.11) can be expressed in terms of a (VG) pro- 
cess with random scale Mg^- specified by (2.9). If for R-r = k{T, Ut), Rt-Y is 
bounded then one can obtain perfect samples of the distribution of T(r), 
via [15]. For certain k, that are not necessarily bounded, one can use the 
density formula in (iii) of Theorem 2.2. 

The next result introduces a nonrandom time change that leads to a 
significant reduction in complexity. First, define for m > 0, 

1 



4>m{lj) = \/2 + fi'^m'^/m + fi and 6rn,(/^) 
and Cm(/u) = fem(/^)/<Am(/^)- 



my^2 + /x^m^ 



Theorem 2.3. For the convoluted subordinator in Theorem 2.2, the 
time changed process {Xq{s) : s > 0) := (W^(T((1 — e~*)/0)) : s > 0) satisfies 
for each fixed s > 0, 

where Mi^s* satisfies (2.10) for t = s* = {1 — e~'^)/9 andp=l—e~'^. 

(i) Furthermore, for each fixed s, Xg{.s) given Mi^g^ = m? , for m > 0, 
follows a double exponential distribution, with density 

f I I N / ^m(/u)e^'^'"(^\ z<^, 
Ixois) [^m-<^ 6™(;u)e-^^-(-^) , z>0, 

and distribution function 

TP I I N _ /cm(/i)e^['^'"(^)], z<0, 
'^^"'"'^ - 1 cM + CrX-m - e-^[^-(-'^W), z > 0. 
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(ii) Hence, if the price process in (2.11) is based on substituting T[t) with 
the time time changed clock, T{{\ — e~^)/9) for s > 0, then the fair price is 
given by 



Substituting t = (1 - e~*)/0 = s*, p = 9s*, yields a GGC(p, Rs*Y) variable, 
which is also a GGC {1, Rs^YS^p). Statement (i) is straightforward. Statement 
(ii) is also not difficult to verify. □ 

In order to evaluate the price in (2.12), it remains to evaluate Mi^^* ■ We 
sketch out the details to do this via the Double CFTP perfect sampler in 
[15]. The deterministic time change allows us to exploit generally the most 
efficient case, 61 = 1, of the Double CFTP. 

First, note again that, R-r* = k{T*,UT*) where k{t, y) is a known function, 
and [/ is a Uniform[0, 1] random variable. Hence, in order to sample Rr* we 
simply need to draw U. Note that we write U to distinguish it from the 
uniform variables we introduce below denoted as C/j. Assuming 



is bounded by a positive constant c, the Double CFTP exact sampler in [15] 
is based on the following steps: 

Backward phase. For i = —1, —2, . . .: keep generating {Ui, Di, D'^) and stor- 
ing (A, A') until Ut < \Dj - D'j\/{2c). Keep T. 
Set starting point. Set Mi^^* = Df A Dj + 2cC/t. 

Forward phase. For i = T + 1, T + 2, . . . , — 1: given (Ai A', Mi^,-*) previ- 
ously stored, do the following step: generate U' uniform [0, 1], ^i/2) a-iid gen- 
erate U uniform [0, 1], and construct X={1- U)Mi^r* + f/A6/2 + UD[{1 - 
ii/2)- Repeat this step until: 



where, r* 



(2.12) 



e"''^E[(5^ - A")+] =E[DE((T2Mi,^*,i^,r)], 
(1 - e~^)/e, and for z = \og{So/K) + rr, 
DE(y2,K,r) = 5oAi/2(2|y) - e-^-^KA/2(z|y). 



Proof 



The result follows from the fact that, again, 
^T(t) (w) = tE[V^L,(a;fc(t, Ut))] = 0tE[log(l + u:RtY)]. 



D = Rr'Yip = k{T*,UT*)YCj 




1 



Ml,,. 



> 1/c 



Output. Return Mi,,-*. 
See [15] for more details. 
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Remark 2.2. These results, which are considerably simplified by using 
the deterministic time change, apply to a wide choice of kernels. It would 
also be nice to find models for T{t) whose marginal distributions were not 
strongly dependent on the form of the kernel. Even better, would be the 
ability to specify laws in a manner similar to how one selects the BDLP 
of an OU to induce general self-decomposable laws for the instantaneous 
volatility v{t). In the next few sections, we will show that this can be done 
for convoluted subordinators we refer to as quantile clocks. 

3. Quantile clocks. As in the Introduction, let R denote a positive ran- 
dom variable with continuous strictly increasing cumulative distribution 
function Fr. Let Qr denote its corresponding quantile function, that is, the 
continuous inverse of the cumulative distribution function. Furthermore, as- 
sume that E[i?] < oo. Then for a subordinator L, we say that the process 
Tr = {Tnit) : t > 0), defined as 

(3.1) TR{t):= J^Qr(^(^1-^^ ^L{ds) for t > 

is a quantile clock with parameters {R, L). Note here that R does not depend 
on t. Furthermore, Qr can be evaluated numerically in many cases, even 
though it may not have a closed form. 

Proposition 3.1. A quantile clock Tr = {TR(t) : t > 0) with parameters 
(R, L) has the following properties: 

(i) If the support of the density of R, say fR, is of the form [0, 6), 6 > 0, 
that is, (5i?(0) = 0, then Tr are random processes with samples paths that 
are almost sure strictly continuous and strictly increasing. 

(ii) Suppose the density of R has support starting at a > 0, that is, 
Qr(0) = a, then there is a positive random variable R with Q^(0) = 0, such 

that R = R + a and Qr{u) = Qp^{u) + a for u G [0, 1] . Hence, it follows that 
the clock can be represented as 

(3.2) TR{t) = T^{t) + aL{t), t>0, 

where satisfies (i). Note Tp, is an {R,L) quantile clock and is obviously 
not independent of L. 

(iii) For each fixed t, the marginal distribution 

TR{t)^at), 

where C, is a subordinator such that C,{1) is a random variable with Laplace 
exponent 



V'c(u;)=E[V'L(u;i?)] = VTH(i)(^)- 
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(iv) That is, the Levy density of C, has the form 

/•oo 

Pd^)= PLis/r)r^^Fii{dr). 
Jo 

Note furthermore that for a constant c, TcR{t) = cC,{t). 

Proof. Setting k{t,s) = Qr{{1 - f)+)) it follows from [4], Proposition 
1, that in order to verify statement (i) we only need to check whether k{t, s) 
satisfies conditions (a), (b), (c). Conditions (b) and (c) are obvious and it 
remains to check the integrability condition, which follows from 

j^Qni^-j^ ^ds = tj^ QR{u)du = tE[R]<oo, 

since Qji{U) = R. Now using this, and standard results for linear functionals 
of Levy processes, we see that for each fixed t, the Laplace exponent of Tii{t) 
is given by 

^ V'L^wQiJ^^l-l) y^ds = tj^ ^L{i^QR{u))du = t¥.[il:L{ujR)\ 

verifying (iii). Statements (ii) and (iv) follows easily from (i) and (iii). Note 
that the quantile function of R in statement (ii) violates condition (b). □ 

We now give an interesting example that has explicit laws. 

Example 3.1 (Arcsine/Bessel occupation time quantile clocks driven by 
L = 7e). First, recall the exponentially tilted stable subordinator Sa dis- 
cussed in Remark 1.2. Suppose that one specifies R = /3i/2,i/2 and L = 751, a 
gamma(^) subordinator. Then 

<3(/3i/2,i/2)N = sm^(|n^, 0<w<l, 

and the quantile clock is defined as 

for t > 0. It fohows (see [12, 20, 22]) that for each fixed t, 

^/9i/2,i/2(*) = -5^1/2 (726»t/2) = letl^et+i/2,et+i/2- 

More generally, for each fixed < a < 1, let Oq(s) = 1{^Bu>o) du denote the 
time spent positive up till time s of a symmetrized Bessel process {Bu,u > 0) 
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of dimension 2 — 2q!, see [1]. Then setting R = 0^(1) := the quantile of 
0„ is 

l/a 



Qoc (u) = '^-'^"j^^l where Qx^ (u) 



Sm TTQU 



sin(7ra(l — u)) 



is the quantile function of the ratio of i.i.d. positive stable random variables 
Xa = Sa/S'^. Then, from James [20] (see Section 7), the clock T©^ with 
parameters {Oa,Je), satisfies for each fixed t, 



where 



Oa,et = Pet,iOa,et + {l-Pet,i)Oa 

= /3a+et,l-aOa,a+et + (1 " /^a+et,l~a)Cl/2 

are random variables corresponding to the time spent positive of general- 
ized Bessel bridges as explained in [20], Section 5. These variables can be 
exactly sampled in various ways as explained in [15]. Furthermore, from [20], 
Proposition 5.3, it follows that for < p = 6t < 1, 

To^{p/0) =7lOa,p, 

where Oa,p = /3p,i-pOQ,p is GGC(l,OaCp) with density 

f- M-^.^-isin^^arctanf^i^^l^^^^^^^'l^ 

X [y2" + 2y-(l-yrcos(a7r) + (l-2/f"]-P/(2"), 

< y < 1. In general, the process {W^{TQ^{t)),t > 0) has almost surely con- 
tinuous sample paths and satisfies, for each fixed t, 

E[e-^M(^o.W)] = 2^*/-(l + (1 + (a;V2 - i/ua;))°)~'*/". 



Remark 3.1. The last example shows that the class of quantile clocks 
where L is a GGC{9,Y) subordinator is equivalent in a marginal sense to 
the representation of GGC variables in terms of Wiener-Gamma integrals 
as defined and presented in [22]. That manuscript, along with the works 
of [5, 19-21], yield many examples of quantile functions which can be used 
to construct quantile clocks with explicit laws, of which quite a few are 
constructed from Qxa- We shall encounter some more examples in Section 
6, although in a slightly different context. 
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4. Choosing quantile clocks to have specific laws in Us- The results in 
the previous section suggest that the the marginal distributions of the quan- 
tile clocks, while equating nicely to the marginals of a subordinator ^, are 
strongly dependent on a random variable R, induced by Qr. Noting that 
Qr{u) is in fact a deterministic function one would like to be able to choose 
explicit laws for Tr, regardless of the function Qr. That is to say, how does 
one choose L so that Tji{t) has a marginal distribution not obviously de- 
pending on R7 For example, for each Qr how does one choose L so that 

TR{t) = Sa{t)'^ Or how does one choose Tji{t) so that a log price process 
W^{Tji{t)) has marginal distributions that are equivalent to a CGMY pro- 
cess? Finally, for different quantile functions Qr^, Qr2, Ri not equivalent 
in distribution to R2, how to choose the driving Levy processes, say Li and 

L2, such that marginally for each fixed t, TR^(t) = Tji^{t)? 

We saw that this was difficult in the case of general convoluted subor- 
dinators as their laws depend strongly on t through the kernel or variable 
k{t,Ut). However, Proposition 3.1 shows that one can represent 

(4.1) TR{t)^ C QR{y)L{tdy), 

Jo 

and there is a clear separation of the effects of t and Qr. This is similar 
to the case of the OU models f (0), see [2, 3, 26], where every positive self- 
decomposable random variable can be represented as 

-00 Jo 

where i9, is a subordinator referred to as a OU-BDLP. More strikingly, there 
is a simple way of obtaining any desired self-decomposable law for v{0) by 
choosing the BDLP according to either of the equations 

(4.3) Vi9(w) =a;V'^,(o)(^) and p^{x) = -p^(^o){x) - xp'^^^^^ix). 
We noticed from (4.1) that ii R = U^'^ = /3s^i for 6>0, then 

Jo 

where, we substitute Z for L, and ^^^^ are subordinators having laws in Ju- 
rek's [17, 23-25] Us class of generalized s-selfdecomposbale laws, where Us C 
J'. The case of (5 = 1, corresponds to Jurek's U = Ui class of s-selfdecomposable 
class. Using Jurek [24, 25], one sees that for < (5i < 1 < ^2 < co. 



g+cG CCcUs.CU CUs^Cj. 
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It follows that for each C^^^ G Us there is a, Z such that 



and hence from [14], Lemma 1, which can be verified directly by taking 
derivatives with respect to 00 of 

(4.4) w'^Vc(^) (w) = / il)z[u)5u^-^ du, 

Jo 

one sees that the Z is related to (^^^^ , and hence T^^ ^ , by the 

equation 

(4.5) ipz{(^) = '4^<^(s){u}) + ^uji;'^^s){uj)- 

This is analogous to the relationships between v{0) and its OU-BDLP •d, 
given in (4.2) and (4.3). We shall show that this relationship becomes more 
explicit as one restricts their choices of laws for <^^^^ to £, Q and Q^. 

Remark 4.1. The specifications in (4.5) and its refinements now allow 
us to specify any law in Us for quantile clocks based on Qjji/s, analogous to 
the case of the BNS-OU v(t). This, as far as we know, is the first instance 
where such a property has been noticed for convoluted subordinators. How- 
ever, in terms of choices of Qn this is still restrictive. The next results show 
how, for a large class of quantile functions Qr, to choose L such that for 

each fixed t, TR{t) = C^^\t) G U5. 

Theorem 4.1. Consider the specifications for a quantile clock 



t„W = ^*c3„((i-;)Jl(<!, 



with parameters {R, L) . Now select R so that its density has hounded support 
and let Y denote a positive bounded random variable such that 

(4.6) RY = U^/^ = (3s^i 

for a fixed 5 > 0; if R = cU^/^ , then Y = 1/c. Suppose that one wishes to 
choose L such that 

Tn{t)^C^'Ht), 
where C^^^ is a subordinator with 

^^;^,s,{co)=E[^l^z{u^U'^% 
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where Z is aUs-BDLP, satisfying (4-5), and hence C,^^^ is in Us- Then for 
this Z, L is chosen such that 

(4.7) VlI'^) =IE[^z(wy)] equivalently pl{x) =¥.[pz{x/Y)Y-\ 
That is, 

(4.8) ^l(^) =E[^^(,)(a;y)] + i^E[yv4w (^>^)]. 
Note that Y is chosen independent of R and Z. 

Proof. The difficulty of tliis result is envisioning its construction. The 
proof itself is otherwise straightforward, since (4.7) and (4.5) implies that 

i:^(,){u:)=¥.[,ijL{ooR)]=¥.[iljz{ooU^'% □ 
We now specialize this result to self-decomposable laws. 

Theorem 4.2. Consider quantile clocks Tpi with parameters {R,L) sat- 
isfying (4-6) and (^^.7), (4-8)- The next result describes further specifications 
in order for Tr,{(^^^) to have laws in C,Q and respectively. 

I. Tr G C: IfTji is selected such that its marginal laws are self- decomposable, 
then it is known that there exists a subordinator such that 

e-^(y-'^i^{Xds)= / e-'^{ds)=v{0). 

-oo Jo 

Furthermore, adapting (4-3), one has 

(4.9) ^/Ji){uJ) =uj'il;'^^s){^) and p^{x) = -p^(s){x) - xp'^^s){x). 

(i) Hence, the BDLP L has to be chosen such that the Levy density of 
Z is 



= p^(6){x) + ^p{>{x). 

That is, 

Pl{x) = E (^p^,s, (x/Y) + ^p4^/Y)^ 

(ii) Statement (i) implies that the subordinators are related as follows: 
Z{t) = C^^\t)+^{t/5), t>0. 
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II. Tr G Q : IfTji, (C^^^) is selected such that its marginal laws are GGC{u), 
it follows that 

/•OO /"OO 

p^(s){x) = x~^ / e~^^u{dy) and p^{x) = I e~^^yi'{dy). 
Jo Jo 

(i) Hence, the Levy density of Z, say pz, satisfies 

/•OO 

pz{x)= / e-^y[x-^ + y/5]u[dy). 
Jo 

(ii) Equivalently, 



log(l+a;/y)+ ^ 



5 y + uj 



v{dy). 



III. Tr G Q^: IfT^, iC^^^) is selected such that its marginal laws are GGC{6, V), 

(i) then L must be selected such that it is equivalent in distribution to 
the subordinator 

L{s) = Cs,Yis) + ^Y{s9/S), s>0, 

where Cs,Y is a GGC{9,VY) subordinator and 

N[s) 
k=l 

where are independent exponential{l) variables, {Y^) are i.i.d. 

variables with distribution Fy , (Vfc) are i.i.d. Fy , and N(s) denotes 
a homogeneous Poisson process with E[A^(s)] = s. 

(ii) As special cases Tji{t) = 7e(t) is obtained by setting V = 1. 



The proof of this result is fairly immediate from the definitions of the 
various classes, details are omitted. 

Remark 4.2. Theorem 4.2 shows that in order to specify Tr to have 
laws in C, one only needs to identify the OU-BDLP i? that leads to a corre- 
sponding stationary law for v{t) =Tr(1) = and use it appropriately 
to define L. One may consult for instance [3] for many explicit examples 'd, 
and the laws they induce. 



5. Choosing R and Y such that RY = U^^^ . The results in the pre- 
vious section show that for a deterministic quantile function Qji one can 
choose quite arbitrary marginal laws for Tr, analogous to the case of f (0), 
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provided that one identifies a variable Y such that RY = U^/^ = fi^ i for a 
fixed 5 > 0. Notice that 

(5.1) Q^,i,{u) = [QR{u)f'. 

The easiest case is to choose R = U and 1" = 1, which as seen from (5.1) 
leads to quantile clocks corresponding to the Holmgren-Liouville convoluted 
subordinators discussed in [4]. The equation (5.1) suggests that one may 
always work with the pair satisfying the solution RY = U and then obviously 
^i/5yi/5 ^ ^1/5^ However, the case of 5 = 1 may not always be the most 
obvious. 

Example 5.1 (Beta variables including the arcsine distribution). Con- 
sider the case of products of independent beta variables 

(5.2) (3s,..sf3.,i+s~. = f3s,i = U^" 

for 5 < K <1 + 5. Hence, for each (5, one can choose many [R,Y), ranging 
over (5 < K < 1 + 5, such that 

{R^'\Y^I')^{Ps,.-5,(in,i+5-.) or 

(5.3) 

(i?l/^yl/^) = (/3,,l+5„«,/35,K_5). 
Furthermore, for some 6 > 0, and each fixed 5 the variables in (5.3) lead to 
variables R^/^ and y^/^, not having beta distributions, that satisfy (i?y)^/'' = 

Lets look at a special case of this in more detail. 

Example 5.2 (Kumaraswamy and generalized arcsine clocks). Setting 
K=l and 5 = a, (5.3) leads to the choice of the pair 

(5.4) (/3,,i_,, 1 - C/i/-) = (/3„,i-a,/3i,a), 

such that i?i/"yi/'^ = [/^/" = /3q,^i, where the first component in (5.4) has 
the generalized arcsine law which arises in many studies of random processes. 
Setting R^/^ = [1 — U'^l°-^^^ = Ka,b leads to the quantile function 

of a Kumaraswamy distribution. Hence, the law of a Kumaraswamy quantile 
clock Tk^^, that is, with parameters {Kafi,L), can be specified such that 
its marginals satisfy 
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where is a subordinator having any law in lAab and hence in C Specif- 
ically, this is done by the choice of 

V'l(^) =E[V<^(.,)(w(/3,,i„,)^/'')] + ^wE[Vj(.,)(a;(/3,,i_,)^/^)(/3,,i_„)^/']. 

Note that if we instead choose R = /3i/2,i/2 and hence R^/^ = (/3i/2,i/2)^^* 
we obtain clocks based on the arcsine law with quantiles 

(5.5) Q^,/,(n):=[Q(^^^^^^^)(^)]V^ = sin2/''(^n). 

This case can be compared with Example 3.1. 

5.1. Selections based on decompositions of an expoiie'n.t[al{l) variable. In 
general, by rescaling to [0, 1] , we see that choosing an R and Y to satisfy 

(4.6) for some (5 > is equivalent to choosing variables iji and iy such that 

(5.6) ij^ + iy£j^/S. 

Recall also the relationship between quantiles of a positive variable X and 

Q^_;,(n)=e-Q^(^-"\ 0<u<l. 

There are obviously many pairs satisfying (5.6). We next look at two 
different types of examples based on suggestions made to us by Prof. Marc 
Yor. 

Example 5.3 [Fractional and integer parts of an exponential(l)]. We 
first note that one of the reasons the following example is interesting is that 
it identifies a concrete example of a quantile clock that is of the form 

TR{t) = T^ + aL{t) 

as specified in statement (ii) of Proposition 3.1, but where we can apply 
Theorem 4.1. Now, following Chaumont and Yor ([11], page 42, Exercise 
2.18), let [71] and {71} denote the fractional part and integer part of an 
exponential(l) variable 71, then (remarkably) these variables are indepen- 
dent and obviously satisfy 

[71] + {71} = 71- 

In this case, RY = U, for 

(5.7) R£e-^^^^ = U{l-e-^) + e-^ and Y = e-^^'\ 

where {71} is a geometric random variable with success probability 1 — e~^ 
and values in {0, 1,2,...}. We say such a variable is geometric (1 — e~^). We 
can extend this case as follows. 
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Proposition 5.1. For 0<p<l, letUp = Up+ {l-p) and let Xp be a 
geometric(p) variable. Then 

Proof. It is easy to verify that the Laplace transforms of —\og{Up + 
(1 — p)) and Xp[— log(l — p)] are given, respectively, by 

and 



p(l + a;) 1 - (1 

Hence, their product is 1/{\+uj), which is the desired result. □ 



That is, there are variables UpYp = U for 



R = Up = Up+{l-p) and Yp = e 



d 



for Xp a geometric(p) variable for each < p < 1. Naturally Ui = U . Hence, 
for each fixed p, 

T^^ {t) = L{t){l-p)+p ^L{ds) for t > 0. 

We now state a result which now follows obviously from Theorem 4.1 and 
applies for quantiles clocks based on the variable tjp^^ . 

Proposition 5.2. For each 5 >Q, and < p < 1, set \ = — log(l — p), 
then the quantile clock 



T^i/s (t) 







{l-p)+p{ 1-J 



1/5 



L{ds) 



can be specified such that for each t, T-i/s{t) = C^^^i^) 'if L is chosen 

Up 

such that 

where Xp is a geometric(p) random variable. Whenp = 1, the quantile clock 
is continuous, otherwise it has jumps. 

Example 5.4 (Splitting the Laplace exponent of 71, part I). Next, con- 
sider the Laplace exponent of 'yi/5, 

/•oo 

log(l+w/(^)= / {l-e-"^/^)s-h-'ds, 
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then choose £r and iy according to the decomposition of the Levy density 

s"^e"* = 7ri(s) +7r2(s). 

In particular, iji and £y are infinitely divisible based on the Levy densities 
vTi and TT2, respectively. The simplest case is where 

s'^e"** = (1 - a)s~^e~'' + as'^e"** 

leading to 

(^y)l/5^[g-7.e-7i-.]V5 d f^l/<5^ 
So, for instance, the quantile clock with parameters (e"^",^), that is, 

Te-rcit) = re-Q^"(i-(i'^/*)+)L((is) 
Jo 

is based on a nontrivial quantile which can however be evaluated by vari- 
ous computational packages. Furthermore, our results show that despite the 
complexity of this clock we can choose quite general marginal laws for T^-ja , 
not directly depending on the quantile function, by working with a BDLP 
satisfying 







r(l-a) 



6. GGC decompositions of a 7 subordinator and resulting clocks. We 

now identify a very large class of variables satisfying (5.6) based on Example 
5.4, using variables in B and in particular, the GGC class 

Theorem 6.1. Let 0.^ he a subordinator in B with Levy density pugis) = 
e~*^g(x) dx, where q{x) is a nonnegative measure such that q{x) < 1 for 
x>5. Then, there is another subordinator in B, say Vts, with Levy density 
Pn ~ Is ^~^^[^ ~ Q{x)]dx, such that the Levy density of a gamma(l) 
subordinator with scale 1/6, has the decomposition 

(6.1) s-h-'' = p^^{s) + pnM- 

Hence, the gamma subordinator can be expressed as a sum of the subordina- 
tors o.nd 0^, which implies for each fixed 9 

Proof. The Levy density of the subordinator {'yi{t)/5:t > 0) can be 
expressed as 



oo 



5 
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leading easily to (6.1). □ 

We next describe an interesting special case involving variables in Q^. 

Theorem 6.2. Assume that Us in Theorem 6.1 is a GGC{l,V/6) sub- 

ordinator for V a random variable in [0,1], and let X = {I — V)/V. Let 
(Et(F) :t>0) denote a subordinator with Levy density denoted as PEi(y)- 
Note that T,i{V) is not random in V. 

(i) The Levy density of a gamma(l) subordinator with scale 1/5 has the 
decomposition 

where 

/oo 
e~'''[l-Fi/v{x)]dx 

(6.2) 

= -e-"(l-E[e-^^]). 
s 

(ii) If < Q = E[-log(y)] < oo, then Pt.^(v){s) = QfA4X){s), where 
/Ae(x)(s) is a density of a random variable denoted as Ae{X), determined 
by (6.2). In this case T^tiV) is a compound Poisson process representable as 

N{gt) 
k=l 

(Afe) are i.i.d. random variables equal in distribution to Ae(^). 

(iii) In general, for each fixed 0, 

(6.3) ^e = ^e{V) + leMe, 

where Me = Pe^Me + (1 - /3e,i)y. 

Proof. From [6], Section 9, we know that corresponds to a GGC(1, V) 
subordinator if q{x) = Fi/y{x). Hence, scaling by 6, and using known prop- 
erties of variables in concludes the result. □ 

Note from, for instance, [6], Example 9.2.3, it follows that for each 5 > 0, 

one can choose ^5(1) = —log{/3s^K~5) or Us{l) = — log(/3K_i+5_K) for the beta 
variables in Example 5.1 satisfying (5.2). Furthermore, among these, the 

only choice corresponding to a GGC variable is 0,s{l) = — log(/35^i). How- 
ever, Theorem 6.2 allows us to construct many quantile clocks based on 
variables in whose distributional properties are explicit. We next de- 
scribe an interesting property of the variable Si(p). 



QUANTILE CLOCKS 25 

Proposition 6.1. Suppose that Tn is a quantile clock with parameters 
{R, L) such that, for an independent variable Y , 

RY = U^'^ = f3s,i = UP 
for 6 = 1/p > 1. If L is chosen such that 
(6.4) i;L{i^) = E[ijz{ujYe-^'^P^)], 

where Z is a Ui-BDLP, satisfying (4-5) for a suhordinator ^ £Ui, then 
^i?(i) = C{t) for each fixed t. 

Proof. Setting V = p, it follows from (6.3), with = 1, that 

71 =Pli + ^i{p), 

which gives the identity 

C/Pe-^i(p) = [/. 
Hence, (6.4) leads to E[V'L(wi?)] = IE[V'z(wt/)]. □ 

6.1. Interpreting Titiy) via diffusions straddling an exponential time. Pro- 
vided that < E[— log(y)] < oo, the random variable Ae(-'^), with density 
defined by (6.2) has an interesting interpretation that we now discuss. This 
will also give us an opportunity to describe some more explicit examples of 
Qn- Let e/X denote an independent exponential 1) time e divided by an 
independent variable X with distribution characterized, for bounded mea- 
surable functions H , by 

E[H{X)]=E[H[X)\og{l + X)]/Q. 

Now let {Ti!"s'^'' , s > 0} denote a recurrent linear diffusion starting at whose 
inverse local time, in this case, is a gamma(l) subordinator. Define for any 
t>0, 

gt := sup{s < t; TZf'^^ =0}, dt := inf{s > t, TZf'^^ = 0}. 

Then given X = X,it follows from (6.2) that for an independent exponential(A) 
variable e/A, the random variable 

Ae(A) = (ie/A -C/e/A 

corresponds to the length of excursion of TZ^^'^^ above straddling an 
exponential(A) time. See, for instance, [34], Section 4, for this description for 
more general TZ as well as [5, 32, 36]. In addition, see [5, 19, 22] for Ae(A) 
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representation as a variable in Q^. Hence, Ae{X) interprets as Ae(A) but 
now for a random time e/X , with c.d.f. F^/x satisfying 

l-^e/x(y)=E[e"^^log(l + X)]^. 

It follows that for X = (1 — p)/p and g = — log(p), that 

N{et) 

k=l 

where (^/a'S'I/a) i-^-^- ^^P^^^^ °^ (4/A,fe/A)- 

6.2. Some related examples. From the results in [5] (see also [19, 22]), we 
can consider more generally TZ^'^'^\ in place of TZ^^'^\ which, for < a < 1 is 
now a process whose inverse local time is distributed as a generalized gamma 
subordinator with Levy density specified by s~"~^e~'*/r(l — a) for s > 0. 

Furthermore, for A = 1, A^'^\ is the generalization of Ae(l) = Ag'''^^ with 
density 

(6.5) Ai-^) I ax---^e-^il-e-p ^ 
^ ' [2"-l]r(l-a) 

Note that the variable Ua,e = e has density 

(6.6) fu (,)^«H°g(")r7^(l-^) forO<n<l. 
V ; ju^,^y ) [2"-l]r(l-a) 



= - 1 



In addition, [5] show that Ae ' is GGC(1 — a,BQ,), where Dq, satisfies 

1 



with 



log(Xi„,) = log(5i_,/5i„ J = log(G,/(l - G«)). 

\ — a 



Furthermore, G1/2 = /3i/2,i/25 Gi = [/ and I/Gq = 1 + e'^'' for r/ a standard 

Cauchy variable. Furthermore 7i_qC/ = 7i/3i_Q^i+a is GGC(1 — a,Ga). We 
now look at some special case of Theorem 6.2. 



Proposition 6.2. LeiO<a<l. 
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(i) Then for V = Dq, and X = Ga, 

(6.7) =A^')+7,M„ + Si(B,) 

where A^"'^^ has density (6.5). = /3a,iM„ + (1 - (3a,i)na- When a =1/2, 

71/2-^1/2 = a2'^'^'^\ otherwise, jaMa has an explicit density given in [19], 
Theorem 

(ii) ForV = Qa^ and X = (Xi^^.)^^""^/", 

71 = 7i_Q,C/ + i;i_„(Ga) + 7« 

(6.8) =7i_«C/ + 7„M„ + Si(G,) 

= 7i_a[/ + Si(^i_q,Gq,), 

w/iere = /3a,iMa + (1 - /Sa.O'^a- 

Proof. Note that aI'^'^'' +7oMa is GGC(1, D^) and the second equality 
above is a direct consequence of Theorem 6.1. Additionally, it follows that 

since A^"'^^ is GGC(1 - «,©„), 7i_„ = Ai""'^^ + Si_„(D«)- The last equality 

follows from aI"'^^ is GGC(1, 6-aID)a)- □ 

Recall the occupation time variables in Example 3.1, where the variable 

Oq^p has density (3.3), and additionally Oq = 1/2 and Oi = ^1/2- Then using 
Theorems 6.1 or 6.2, we obtain the following result. 

Proposition 6.3. Let < a < 1 and < p < 1. 

(i) Then for V = Oq, and X = X^ , 

(6.9) 7i=7id„,p + Sp(Oa)+7i-p- 

(ii) When a = l/2, F = /3i/2,i/2 a^^c^ ^ = 71/2/7^/2' 
(6-10) 71 = 7p/3p+i/2,p+i/2 + Sp(/3i/2,i/2) + 7i-p- 

Remark 6.1. The diffusions above belong to a more general family, 
7^(a,b) for < a < 1, 6 > with inverse local time corresponding to a gen- 
eralized gamma subordinator with Levy density Cs~"~^e~**. Hence, the 
density of variables A^"'^^ is proportional to s~"~^e~^*(l — e~^). In par- 
ticular, for b = 0, the variable Ai"'°^ = 7i_c,/C/i/° ^ GGC(1 - a, 1/G«) 
variable. See [5, 19, 22] for more details. 
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Remark 6.2. Reference [5], Theorem 1.4, yields the following decom- 
position of 7i , for < a < 1 : 

71 =7iGa + 7jGi_Q,. 
As a special case, with a = or 1, 

7i = 7if/ + 7;/(l+e"''). 

Combining this fact with with (6.3), with V = Gq = 1/(1 + e'^^), leads to the 
interesting identity 

Si(Go) =71^0 

with Levy density given by (6.2). This follows since 71 f/ is a GGC(l,Go) 
variable. 

7. Composition of quantile clocks. We now highlight an important prop- 
erty of the general class of quantile clocks. Recall from Proposition 3.1 that 
for each fixed t, the marginal distribution of a quantile clock Tr, with pa- 
rameters {R,L), satisfies 

TMi) = C(t), 

where is a subordinator such that ^(1) has Laplace exponent 

V'c(a;) = E[V'L(a;i2)] = VT«(i)H- 

An important operation for Levy processes is the composition of Levy pro- 
cesses, in financial applications this is associated with time changed pro- 
cesses. The fact that quantile clocks behave marginally like a subordinator 
allows us to obtain the following results. 

First, we can discuss the composition of two independent quantile clocks 
Tr^, Trj, with parameters (i2i,Li) and (i?2,L2), respectively, which can be 
written as 

(7.1) Tr, {Tr, it) ) = Qr, ( (1 - ^) Li (ds) 

for t > 0. The apparently complicated random process appearing in (7.1) is 
no longer a quantile clock. However, as the next result shows, its marginals 
are easy to describe. We use the notation o to denote the composition of 
functions, so, for instance, Trj o Tr^ means the operation in (7.1). 

Proposition 7.1. Let Tr., i = l,...,k, denote independent quantile 
clocks such that pointwise TR-{t) = C,i{t) for independent subordinators with 
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corresponding Laplace exponents iI:q{uj) for i = 1, . . . ,k. Then the composi- 
tion Tfc := (Tfc (t) = Tr^ o • • • o Trj. {t):t> 0) is an increasing process such that 
for each fixed t, 

where Qk is a subordinator with Laplace exponent 

If each Tr. is a continuous process, then is a continuous process. 

Proof. It suffices to show this for k = 2. But this is immediate from 
Proposition 3.1, since given Tr^, 

i'Tn^ {Ta^ (t)) {^) = Tr^ (OV'Ci ('^)- ^ 

Of course one can also compose these clocks with subordinators as follows. 
The next result is immediate. 

Proposition 7.2. Let Tr denote a quantile clock that satisfies Tji{t) = 
Ci{t) for some subordinator (^i. Furthermore, let Q2 denote a subordinator 
independent of Tr and Ci • Then for each fixed t 

TR{Ut)) = Qi{Ut)) and C2(rij(i)) = C2(CiW)- 

We now illustrate an important special case. 

Proposition 7.3. For < a < 1 and < /3 < 1, one can use the speci- 
fications in Theorem 4-2 to construct independent clocks Trj and Tr^ such 

that marginally TR^{t) = Sci{t) and TR^{t) = Sp{t), where Sa and Sp are 
independent with Laplace exponents [(1 + a;)" — 1] and [(1 + 00)^ — 1]. Then 
for each fixed t 

TR,{TR,{t))^TR,{TR,{t))^TR,{Sp{t)) ^ S^{TR,{t)) I S^pit), 

that is, for each fixed t the Laplace exponent is t[{l + w)"^ — 1] . Additionally, 
the first two compositions can be specified such that the resulting processes 
are continuous, but the latter compositions always correspond to processes 
with jumps. 
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8. Continuous VG, CGMY, NIG and other price processes. Summa- 
rizing, we have demonstrated that quantile clocks Tr can either be cho- 
sen to have strictly continuous and increasing paths or can be expressed as 
Tj^(t) + aL(t), where is a continuous increasing quantile clock. In general, 
quantile clocks have marginals that are equivalent to those of a subordina- 

tor, ^, for each t, that is, Tf((t) = C(^)- Moreover, we have shown that for a 
large class of quantiles Qji we can choose Tr to have any desired marginal 
law in Us, by choosing a random variable Y and the subordinator L in a 
clearly prescribed fashion. Furthermore, our results in the last section show 
that composition operations involving quantile clocks or quantile clocks with 
subordinators are marginally equivalent in distribution to compositions of 
subordinators. All these properties make them highly desirable components 
in pricing models based on time changes. For example, the processes 

{%{TR{t)):t>0) and {W^{TR,{TR,{t))) -.t > 0) 

can be chosen such that they are processes with continuous trajectories, but 
have simple and familiar marginal laws. In addition, for a subordinator ^, 
the processes 

(W^{aTR{t))):t>0) and {W,,{TR{Cm :t > 0) 

have jumps, exhibit volatility clustering, and otherwise may be chosen to 
have familiar marginal distributions, in fact the same marginal, for many 
choices of Qr. We illustrate these points through some examples that equate 
these processes marginally with some of the most popular Levy processes. 

Example 8.1 (Continuous variance gamma processes). As a first exam- 
ple, it follows from III of Theorem 4.2, that if RY = U, then for each 6 > 0, 
a quantile clock Tj^i/s with parameters {R^^^ ,Ls), can be chosen such that 
for each fixed t, 

%iTRi/sit))^W,i-feit)), 

that is, it has marginal distributions equivalent to the log price of a vari- 
ance gamma (VG) process [28], not depending on 6, if for each 6 > 0, the 
subordinator (depending on 6) 

N{es/s) 

Ls{s) ^ Cs,Y^Ms) + E /^^^fc^'' ^^0' 

k=l 

where C^yV* is a GGC(0,yi/^) subordinator. 

We next show how to obtain price processes whose marginal laws are 
equivalent to a Carr-Geman-Madan-Yor (CGMY) process [7] but otherwise 
possesses continuous sample paths. 
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Example 8.2 (Continuous CGMY processes). For this example, we fol- 
low the exposition in [29]. Let 

, G-M ^ ^ G+M 

A = and n = , 

2 2 ' 

then the Levy density of the log prices of a CGMY process, say xcgmy, is 
given by 

/oi\ / \ r(a)r(l - a) Ax-B\x\i r 

(8-1) Pxcgmy(i)(^)= Vil + a) ^ '"^l for -oo<x<oo 

and < d = 2a < 2. Madan and Yor [29], show that the log price of a CGMY 
process has an explicit representation in terms of a time changed brownian 
motion, XCGMy(*) := ^a(C(0)i where C is a subordinator with Levy density 



.2) 



p^^s) = ^I^e(^'-^')^/2,— iE[e-'^(^V2)(7./7i/2)] 



for 



^.3) V=(aMG + B^^\ /2. 

V 7i/2 / / 



7l/2 

It is evident from (8.2) that Q ^Q. We now give the specifications for a 
quantile clock to have marginals with Levy density (8.2) hence inducing 
price processes that have the marginal distribution of a CGMY process. 
Notice that 

- sp[{s) = ^£^s-"-iE[((l + a) + .y)e-^] 



(8.4) 

Hence, ^(^{uj) = ¥.[ijz{ujU^I^)], for 



(l + «)p^(,) + ^I^,-E[Fe-^]. 



2-r(a)_,^^^^_,^. 



^.5) pz{s) = (1 + a/6)p^{s) + ^^s-^nVe 



Proposition 8.1. Suppose that T^i is a quantile clock with parameters 
{R,L), such that there exists a variable Y satisfying RY = U^l^ for some 
5 > 0. Then for each fixed t, 



WA{TR{t)) ^XCGMYit) 
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specified by (8.1) if the suhordinator L is chosen such that 
p.(.) = (1 + aM)^^^^«— lEle-^/^'^] 

(8.6) 

where V is defined by (8.3), Ca = E[y"] and is the random variable whose 
distribution is proportional to y^Fyidy). When Y = 1, L := Z satisfying 
(8.5). Note also that K[V°'] is finite only if a < 1/2. Hence, 



s-"E[(F/y,)e 



-sV/Yc,] 



is the Levy density of a compound Poisson process only in the case where 
a < 1/2. 

Proof. The result is a special case of Theorem 4.2 and the specifications 
we derived above. In particular, (8.5). □ 

As a specific example with continuous paths, consider again the Ku- 
maraswamy quantile clock with 

TK,,St) = /V - (1 - (1 - s/t)+f^f'L{ds) 
Jo 

for R = Kp^b = (1 - f/i/P)!/'' = ^l^p. Then for each fixed t, 

WA{TK^,St))= XCGMYit), 

if L is selected according to (8.6) with Y = ji^^i_p and 5 = bp. 
If we consider the arcsine clock using (5.5) with 

then for each fixed t, 

Wa{T.i/, (t))^XCGMY(t), 

^^1/2,1/2 

if L is selected according to (8.6) with Y = {1 — U^)^^^ and 6 = b/2. 

The next two cases are from Section 6.2 involving quantile functions that 
can be evaluated numerically. If we consider the clock T i/s based on the 

,e 

variable Ua,e with density (6.6), it follows from Proposition 6.2 that 

Wa{T i/s{t)) = XCGMY{t), 

^ ct,e 
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if L is selected according to (8.6) with 

Y £ g-pl-a(ID)a)+7<:<]/5^ 

If we consider the variables in Remark 6.2 then this leads to a quantile clock 
based on the variable 71 Gq. Hence, 

WA{T^~jrGc/s{t)) = xcGMvit), 
if L is selected according to (8.6) with 

Other examples using (6.9) and (6.10) are based on the pairs 

(^g-7iOQ,p^g-[Sp(Oc,)+7i-p]-) g^y^fj^ (g-7p/3p+i/2,p+i/2^g-Pp(/3i/2,i/2)+7(_p]-)_ 

Finally, if instead one uses the clock 

t 



(l-p)+p(l 



Lids) 



then for each fixed t, 



T^?A(r.-,i/.(t))=XCGMY(t), 



if L is selected according to (8.6) with Y = e~''''p[~^°s(i-p)]/'5 -where again Xp 
is geometric (p). Hence, for each <p < 1, the resulting process has CGMY 
marginals, exhibits volatility clustering, but also has jumps. lfp= 1, then the 
process is continuous and the quantile clock coincides with the Holmgren- 
Liouville clock discussed in [4]. 

It is evident that the specifications (8.6) for L appearing in Proposition 8.1 
can be modified such that Tii{t) has marginals equivalent to a subordinator 
with Levy density 

p^(s) = Cs-"-^E[e-^^] 

for some positive constant C, where y is a much more general random 
variable. That is, L is specified by 

Pl{s) = (1 + a/(5)c„Cs-"-^E[e-''^/^-] 

(8.7) 

+ (c„/5)Cs-"E[(y/y«)e-^^/^"]. 

As a specific example, we next look at the case corresponding to NIG and 
related processes. 
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Example 8.3 (Processes with NIG and related marginals). For this ex- 
ample, let Sa{t) denote any subordinator with Levy density 

Pa{s)= g-°-^e~^ 
i (1 — aj 

and define the Levy process on M by Xa{t) '■= Wn{Sa{t)). It follows that 

Xi/2{i) :=XNlG(i) 

is a normal inverse Gaussian (NIG) process [2]. If denotes a quantile 
clock such that RY = U^/^ , and L is specified according to (8.7), specifi- 
cally using the Levy density pa, with V = 1, then Wfj_{T^\t)) = Xa{t)- In 
particular setting a = 1/2, it follows that for each t, 

%{Tii/'\t))^XmG{t). 

Thus, yielding processes with continuous trajectories but NIG marginals. 
In addition, choosing a and < /3 < 1 such that a/? = 1/2, it follows from 
Proposition 7.3 that 

x/3(r(")(t)) ^iy^(r(")(5;3(t))) ^XNiG(t), 

corresponding to processes with jumps, dependent increments and NIG 
marginal distributions. Note that when Y = 1, corresponding to the quantile 
clock T^l/s, then similar to (8.5), 

N{as/5) 

(8.8) L{t) = S^{il + a/5)t)+ ^ 7)^^^ for t > 0. 

fc=i 

9. Choosing laws for the short memory kernel. We now apply our results 
for quantile clocks to a convoluted subordinator that [4] refer to as a short 
memory kernel. We note that this convoluted subordinator is not a quantile 
clock. 

Theorem 9.1. Let ^ denote a subordinator with self- decomposable laws 
such that the quantile clock with parameters {U, Z), that is, Tjj has marginals 

Tu{t) = C(0 ^ ^- This is achieved by setting the Levy density of Z to be 

(9.1) pz{x) = -xp'(-{x) = pc_{x) + p^{x), 

where "d is the OU-BDLP of u(0) =Tjj{l) = C(l)- Then for Z satisfying 
(9.1), the short memory convoluted subordinator constructed as 

f,{t) = l^^m(^l,^^^^^Z{ds) 
has the following distributional properties: 
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(i) For each fixed t, the Laplace exponent of the r.v. Tir{t), is given by 



J^p^{oJ) + {t-e)oJ^p'^{uJ), t>e. 
(ii) For each fixed t, the marginal distribution of T£{t) is given by 



Ts{t) 



-C{t), t<e, 
C{t) + i}{t-e), t>e. 



Proof. First, notice that in general, for each fixed t, the Levy exponent 
of the random variable Ti;{t) is given by 

ipz \ wmm 1, 



(9.2) 



^f,(i)(^) 







ds. 



So for t < e, (9.2) can be expressed as 

s) 



* f ( (t 
ipz \ wmm 1, — 







£ 



ds = tE 



t'lp(^ I UJ- 



where the last equality follows from (9.1). For t>e, split the interval [0,t] 
into [0, t — e] and (t — e, t] then (9.2) becomes 

(9.3) eE[Vz(wf/)] + (i-e)V'z(w). 

Now use (9.1) to show that (9.3) is equal to 

yielding the result. □ 

Our result now allows one to choose more convenient laws for which 
allows one to easily apply the option pricing formula of [4], as displayed in 
Theorem 2.1, either by exact simulation or FFT methods. We illustrate this 
in the next example. 

Example 9.1 (Short memory convoluted subordinator with NIG related 
marginals). First, it is interesting to recall from [3] that the OU-BDLP, 

leading to v{0) = 5q,(1), as specified in Example 8.3, has Levy density 

a 



Hence, 



r(l - a 

N{as) 



s-"-Ha + sle- 



i?(s) ^ Saias) + 7fJ„ for s > 



k=l 
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Now using (8.8) with 6 = 1, setting the subordinator 

N{as) 
k=l 



Z{s) = + a)s) + 7fJ^ for s > 



leads to the fohowing marginal behavior of the corresponding short-memory 
model, for each fixed t, 



Teit) = { 



N{a(t-e)) 

Mt + ait-e))+ Yl t>e. 

k=l 



Note however that is a continuous process since Z is an infinite activity 
process that satisfies the conditions in [4]. This leads to price processes 
(2.5) with continuous trajectories that have the following marginal behavior. 
W_i/2{^'^'^e{'t)) is equivalent in distribution to 



.2* 



W_y2[cT'-S^{t)j, t<e, 

/ Niait-e)) \ 

W_y2{Ait + ait-e))) + W'_^/Ja'' Y ^i-" ' * > ^• 



k=l 



If one sets a = 1/2, then this reduces to 



.2t 



W.i/2{a'-Sy2{t)), t<e, 

/ N({t-e)/2) 



W^i/2{c^%/2m-e)/2)) + W'_y2[(^' Y ^1/2 )' *>^- 

k=l 



Hence, in this case, for each fixed t <e, the marginal distribution of the price 
process (2.5) follows an NIG distribution, with scale parameters depending 
on t. 
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